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We present a theory of the spin Hall magnetoresistance (SMR) in multilayers made from 
an insulating ferromagnet F, such as yttrium iron garnet (YIG), and a normal metal N with 
spin-orbit interactions, such as platinum (Ft). The SMR is induced by the simultaneous 
action of spin Hall and inverse spin Hall effects and therefore a non-equilibrium proximity 
phenomenon. We compute the SMR in F|N and F|N|F layered systems, treating N by spin- 
diffusion theory with quantum mechanical boundary conditions at the interfaces in terms 
of the spin-mixing conductance. Our results explain the experimentally observed spin Hall 
magnetoresistance in N|F bilayers. For F|N|F spin valves we predict an enhanced SMR 
amplitude when magnetizations are collinear. The SMR and the spin-transfer torques in 
these trilayers can be controlled by the magnetic configuration. 

PACS numbers: 85.75.-d, 73.43. Qt, 72.15.Gd, 72.25.Mk 

I. INTRODUCTION 

Spin currents are a central theme in spintronics since they are intimately associated with the 
manipulation and transport of spins in small structures and devicesjii^ Spin currents can be gener- 
ated by means of the spin Hall effect (SHE) and detected by the inverse spin Hall effect (ISHE)j^ 
Of special interest are multilayers made of normal metals (N) and ferromagnets (F). When an 
electric current flows through N, an SHE spin current flows towards the interfaces, where it can be 
absorbed as a spin-transfer torque (STT) on the ferromagnet. This STT affects the magnetization 
damping^ or even switches the mag netization^iS The ISHE can be used to detect spin currents 
pumped by the magnetization dynamics excited by microwaves^"— or temperature gradients (spin 
Seebeck effect) 

Recently, magnetic insulators have attracted the attention of the spintronics community. Yt- 
trium iron garnets (YIG), a class of ferrimagnetic insulators with a large band gap, are interesting 
because of their very low magnetization damping. Their magnetization can be activated thermally 
to generate the spin Seebeck effect in YIG|Pt bilayers . By means of the SHE, spin waves can be 
electrically excited in YIG via a Pt contact, and, via the ISHE, subsequently detected electrically 
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in another Pt contactji* Spin transport at an N|F interface is governed by the complex spin-mixing 
conductance G-fiM- The prediction of a large real part of G-fi for interfaces of YIG with simple 
metals by first principles calculations^'' has been confirmed by experiments.-l& 

Magnetoresistance (MR) is the property of a material to change the value of its electrical 
resistance under an external magnetic field. In normal metals its origin is the Lorentz force^^^ 
The dependence of the resistance on the angle between current and magnetization in metallic 
ferromagnets is called anisotropic magnetoresistance (AMR). The transverse component of the 
AMR is also called the planar Hall effect (PHE), i.e. the transverse (Hall) voltage found in 
ferromagnets when the magnetization is rotated in the plane of the film.— i2i Both effects are 
symmetric with respect to magnetization reversal, which distinguishes them from the anomalous 
Hall effect (AHE) for magnetizations normal to the film, which changes sign under magnetization 
reversalj^^ The physical origin of AMR, PHE, and AHE is the spin-orbit interaction, in contrast 
to the giant magnetoresistance (GMR), which reflects the change in resistance that accompanies 
the magnetic field-induced magnetic configuration in magnetic multilayers;^ 

Here we propose a theory for a recently discovered magnetoresistance effect in Pt|YIG bilayer 
systemsJ^i2ii2^ This MR is remarkable since YIG is a very good electric insulator such that a 
charge current can only fiow in Pt. We explain this unusual magnetoresistance not in terms of an 
equilibrium static magnetic proximity polarization in Pt,^^ but rather in terms of a non-equilibrium 
proximity effect caused by the simultaneous action of the SHE and ISHE and therefore call it spin 
Hall magnetoresistance (SMR). This effect scales like the square of the spin Hall angle and is 
modulated by the magnetization direction in YIG via the spin-transfer at the N|F interface. Our 
explanation is similar to the Hanle effect-induced magnetoresistance in the two-dimensional electron 
gas proposed by Dyakonov^^fi Here we present the details of our theory, which is based on the spin- 
diffusion approximation in the N layer in the presence of spin-orbit interactions^^'' and quantum 
mechanical boundary conditions at the interface in terms of the spin- mixing conductance We 
also address F|N|F spin valves with electric currents applied parallel to the interface(s) with the 
additional degree of freedom of the relative angle between the two magnetizations directions. 

This paper is organized as follows. We present the model, i.e. spin-diffusion with proper 
boundary conditions in Sec. [Hi In Sec. IIIH we consider an N|F bilayer as shown in Fig. [1] (a). We 
obtain spin accumulation, spin currents and finally the measured charge currents that are compared 
with the experimental SMR. We also find and discuss that the imaginary part of the spin-mixing 
conductance generates an AHE. F|N|F (Fig. [T](b)) spin valves are investigated in Sec. llVi which 
show an enhanced SMR for spacers thinner than the spin-flip diffusion length. We summarize the 
results and give conclusions in Sec. El 



II. TRANSPORT THEORY IN METALS IN CONTACT WITH A MAGNETIC 

INSULATOR 

The spin current density in the non-relativistic limit 

V = en{v0ff + a(E)v) /2= (Jsxjsyjsz) ={jsJsJl) (1) 

is a second-order tensor (in units of the charge current density jc = en {v}), where e = |e| is the 
electron charge, n is the density of the electrons, v is the velocity operator, a is the vector of Pauli 
spin matrices, and (• • • ) denotes an expectation value. The row vectors jsi = en {vcTi + cTiv) /2 
in Eq. ([1]) are the spin current densities polarized in the ^-direction, while the column vectors 
ji = en {vjff + avj) /2 denote the spin current densities with polarization a flowing in the j- 
direction. Ohm's law for metals with spin-orbit interactions can be summarized by the relation 
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between thermodynamic driving forces and currents that reflects Onsager's reciprocity by the 
symmetry of the response matrix:— 
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where //^ = {^sx, fJ-sy, fJ-sz) — /^oi is the spin accumulation, i.e. the spin-dependent chemical 
potential relative to the charge chemical potential fj,o, a is the electric conductivity, 0sH is the 
spin Hall angle, and "x" denotes the vector cross product operating on the gradients of the spin- 
dependent chemical potentials. The spin Hall effect is represented by the lower non-diagonal 
elements that generate the spin currents in the presence of an applied electric field, in the following 
chosen to be in the x-direction E = E^x = —xdxfJ-o/e. The inverse spin Hall effect is governed by 
elements above the diagonal that connect the gradients of the spin accumulations to the charge 
current density. 

The spin accumulation fig is obtained from the spin-diffusion equation in the normal metal 



(3) 



where the spin-diffusion length A = ^/Dt^ is expressed in terms of the charge diffusion constant 
D and spin- flip relaxation time Tgf For films with thickness dj^j m the z-direction 



(4) 



where the constant column vectors A and B are determined by the boundary conditions at the 
interfaces. 

According to Eq. ([2]), the spin current in N consists of diffusion and spin Hall drift contributions. 
Since we are considering a system homogeneous in the x-y plane, we focus on the spin current 
density flowing in the z-direction 

a 



dzfls - JsOV 



2e 



(5) 



4 



where j'^g^ = OsnaE^ is the bare spin Hall current, i.e., the spin current generated directly by the 
SHE. 

The boundary conditions require that jg{z) is continuous at the interfaces z = djsi and 2; = 0. 
The spin current at a vacuum (V) interface vanishes, = 0. The spin current density ji^^ at a 
magnetic interface is governed by the spin accumulation and spin-mixing conductance:^^ 

ejp^ (m) = Grrh x {m x fls) + Gi {rh x flg) , (6) 

where rh = {mx,rny,mz)^ is a unit vector along the magnetization and G^i = Gr+iGi the complex 
spin-mixing interface conductance per unit area. The imaginary part Gi can be interpreted as an 
effective exchange field acting on the spin accumulation. A positive current in Eq. ([6]) corresponds to 
up-spins flowing from F towards N. Since F is an insulator, this spin current density is proportional 
to the spin-transfer acting on the ferromagnet 

Tstt = -^m X (mx jf = ^Jf ) (7) 

With these boundary conditions we determine the coefficients A and B, which leads to the spin 
accumulation 



2eA 

/^s = 

a 



jTv + JsidN)) cosh I + (jfo^y + jP {rh)) cosh 

■(F) 



/sinh^, (8) 



where jK^at) = for F(m) |N|V bilayers and j|(dAr) = —js {rh') for F(m) |N|F(r7i') spin valves 



III. N|r BILAYERS 



In the bilayer the spin accumulation ^ is 

^ ^ ^ ^ osinh ^,p^ , . , 2eA cosh ^ 

nJz) = -yiIi 1^ + jf) {rh) ^ , (9) 

^ ^ ' sinhff- ^ ^ a sinh^ ^ ' 

where //^ = |//s(0)| = (2eA/cj)jfo^ tanh [djv/ (2A)] is the spin accumulation at the interface in the 
absence of spin-transfer, i.e., when G^\^ = 0. 

Using Eq. the spin accumulation at 2; = becomes 



With 



^s(O) = y/iO + — {G, [m (m • /x,(0)) - ^,(0)] + G^m x ;U,(0)} coth ^. (10) 
a A 



rh ■ fls{0) = ruyn'^g, (11) 



^arnxy + rhmy2XGi coth ^ 2AGi coth ^ , ^ 

m X /i, = /i° ^ ^ ^ - /iJO ^ (12 

^ ^ " cj + 2AG, coth a + 2AG,coth^^ ^ ^ 



^ . ^ . n 1 + ^G^coth^^ 

/^.(o) = wh ; ,2 , ' —2 

l + ^G^coth^j +(^GiCoth^^j 

^Gr coth ^ (1 + ^Gr coth + (^Gi coth % 

+ mmy^J.°g- 

l + ^G^coth^^f +^GiCoth^^ 



X n ^G,coth% , , 
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1 + ^G,. coth j + f^Gi coth ^^j 
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the spin current through the F|N interface then reads 

jf ) = t^±rh X (m X y)c7Re '-^ ^ + ^ (m x y)c7lm '-^ ^. (14) 

e ' (T + 2AGt^coth£^ e a + 2AGt^ coth ^ 

The spin accumulation 



pJz) ^sinh^2^ ^„ ^„ ^ 2AGti cosh^ 

„ = —y h m X (mxy)Re + (mx y) Im -f^ 

Ijfis sinh^ L V K y) J ^ + 2AGti coth% sinh^ 



-^s simi 2^ CT -|- zalt|^ cum snm 

then leads to the distributed spin current in N 

f(z) . cosh -cosh ^ ^ 2AGtitanh4f sinh ^ 

= y — - [m X (m X y)Re + (m X y) Im ^ r^- 

jfo cosh|^ L V K yj J ^ + 2AGt^coth^^ sinh^^f 

(16) 

The ISHE drives a charge current in the x-y plane by the diffusion spin current component flowing 
along the z-direction. The total longitudinal (along x) and transverse or Hall (along y) charge 
currents become 



M = 1 + ^SH 

Jc 



cosh^^^ , 2AGtitanh4f sinh 

1^ + (1 - ml) Re ^ 

cosh^ ^ cj + 2AGt; coth sinh^ 



(17) 



JctransC^) .2 , p r , 2AGt4 tanh |f sinh^ 

^ = 0SH {m^ruy Re -m, Im) ^ , — ' (1^) 



K a + 2AGt4 coth sinh 



where jj? = ciit^^: is the charge current driven by the external electric field. 

The charge current vector is the observable in the experiment that is usually expressed in terms 
of the longitudinal and transverse (Hall) resistivities. Averaging the electric currents over the film 
thickness z and expanding the longitudinal resistivity governed by the current in the (x-) direction 
of the applied field to leading order in 9g^, we obtain 

Piong = ^i^ig = (^) « P + Apo + Ap, (1 - ml) , (19) 

C^trans Jc,trans/-^x a i a rnn\ 

Ptrans = 2 — ~ 2 = Apim^my + Ap2mz, (20) 



0", (T^ 
long 



where 



Apo a2 2A dN 

^ = -^SH^tanh-, (21) 

^^^in^Re ^^^ntanh^t ^ (^2) 
P dN a + 2AGt; coth 

^=_,|„Ai„^^^lil!!=!%, (23) 

P div (7 + 2AGt; coth ^ 

where p = o""^ is the intrinsic electric resistivity of the bulk normal metal. ApQ < seems to imply 
that the resistivity is reduced by the spin-orbit interaction. However, this is an effect of the order 
of ^gjj that becomes relevant only when dj^ is sufficiently small. The spin-orbit interaction also 
generates spin-flip scattering that increases the resistance to leading order according to Matthiesen's 
rule. We see that Api (caused mainly by Gr) contributes to the SMR, while Ap2 (caused mainly 
by Gi) contributes only when there is a magnetization component normal to the plane (AHE), as 
discussed below. 
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FIG. 2: (Color online). Normalized ^sx, fJ-sy, jsx, and jsy as functions of z for magnetizations (a) rh — y, 
(b) rh = {x + y) /V2, and (c) m — x for a sample with ~ 12 nm. We adopt the transport parameters 
p = 8.6 X 10^^ rim, A = 1.5 nm, and Gr = 5 x lO^^O^^m^^. For magnetizations rh = y and rh = x, both 
Psx and jsx are 0. 



A. Limit of d = Im G^j, < Re G^j, = Gr 

According to first principles calculations,^'' \Gi\ is at least one order of magnitude smaller than 
Gr for YIG, so Gj = appears to be a good first approximation. In this limit, we plot normalized 
components of spin accumulation (fisx and fisy) and spin current {jsx = Js ' ^ and jsy = jg • v) 
as functions of z for different magnetizations in Fig [2j When the magnetization of F is along 
y, the spin current at the N|F interface {z = 0) vanishes just as for the vacuum interface. By 
rotating the magnetization from y to x, the spin current at the N|F interface and the torque on 
the magnetization is activated, while the spin accumulation is dissipated correspondingly. We note 
that the x-components of both spin accumulation and spin current vanish when the magnetization 
is along x and y, and reach a maximum value at {x + y) / \/2. 

For Gj = the observable transport properties reduce to 

Piong ~ p + Apo + Api (1 - ml) , (24) 
Ptrans ~ Apim^rriy, (25) 
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where 



-- = -^SH^tanh-, (26) 
Api _ ^2 ^ 2AGr tanh^ 

A 



SH^ ^- (27) 



Equations (j24ti25p fully explain the magnetization dependence of SMR in Ref. bsl . while Eq. (j27p 
shows that an SMR exists only when the spin-mixing conductance does not vanish. Since results 
do not depend on the z-component of magnetization, the AHE vanishes in our model when Gj = 0. 

B. Gr > o-/ (2A) 

Here we discuss the limit in which the spin current transverse to rh is completely absorbed as 
an STT without reflection. This ideal situation is actually not so far from reality for the recently 
found large Gr between YIG and noble metals/i^^i^ The spin current at the interface is then 

js c?r>_/( ^ (^jf^ y-j tanh tanh — ^, (28) 

and the maximum magnetoresistance for the bilayer is 

^ = eiH:^tanh%tanh2^. (29) 
p Oat a 2A 

In Sec. IIIIEl we test this limit with available parameters from experiments. 

C. X/dN > 1 

When the spin-diffusion length is much larger than the thickness of N 

/2s(z) A/djv>i ^ ,2z - dN 

— T— = mx{mxy)-y — , 

Ms dN 

while spin current and magnetoresistance vanish. We can interpret this as multiple scattering of the 
spin current at the interfaces; the ISHE has both positive and negative charge current contributions 
that cancel each other. 



D. Spin Hall AHE 

Recent measurements in YIG|Pt display a small AHE- like signal on top of the ordinary Hall 
effect, i.e. a transverse voltage when the magnetization is normal to the film.^*^ As mentioned 
above, an imaginary part of the spin-mixing conductance G, can cause a spin Hall AHE (SHAHE). 

The component of the spin accumulation iigx 

l^s.{z) 2Acosh^ aGn 

n — = -r — ma;m„Re— m^Im — ^ (30) 

a sinh^^f ^ " V + 2AGt;coth^^ ^ ' 
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FIG. 3: (Color online) Calculated Api/p as a function of A for different spin Hall angles 6'sh with (a) Gr = 1 x 
10" ni-2, (b) Gr=5x 10" ni-2, (c) = 10 x lO" m'^, and (d) the ideal limit Gr > cr/{2X). 
The Pt layers are 12-nni-thick with resistivity 8.6 x 10~^r2ni (Sample 1, solid curve) and 7-nm-thick with 
resistivity 4.1 x lO^^'fim (Sample 2, dashed curve). Experimental results are shown as horizontal lines for 
comparisoui^ 



contains a contribution that scales with and contributes a charge current in the transverse (y-) 
direction 

jft"at™^(-2) 9 sinh^ Gtitanh^f , , 

Jc° sinh^ a + 2AG|; coth ^^f ^ 

The transverse resistivity due to this current is 



.(SHAKE) ,1-, 

(SHAKE) Jc,trans / . , 

Ptrans ~ ^ = -^P2m„ (32) 

where 

Ap2 2A26'§h fjGitanh^l^ 2A20|h ciGi tanh^ 



(^(j + 2AG^ coth + (^2AGi coth (^a + 2AG^ coth ^ 



2 ■ 



E. Comparison with experiments 



There are controversies about the values of the material parameters relevant for our theory, i.e. 
the spin-mixing conductance Gf^ of the N|F interface, as well as spin- flip diffusion length A and 
spin Hall angle ^SH in the normal metal. 
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Experimentally, Burrows et al.^^ found for an Au|YIG interface with Go = e^/h. 

= 5.2 X 10^^ m-2; G^"P = 2 X lO^^U'^ m'^ . (33) 

Go 

On the theory side, the spin-mixing conductance from scattering theory for an insulator reads^^ 

^ = iVsh - E <t^ni = Ns^-Y, e<^"^-^"t) , (34) 

" n n 

where ^'ntd) ~ e*''"^*^) is the reflection coefficient of an electron in the quantum channel n on a 
unit area at the N|F interface with unit modulus and phase (5nt(4.) majority (minority) spin, 

and A^sh is the number of transport channels (per unit area) at the Fermi energy, i.e. Ngy^ is the 
Sharvin conductance (for one spin). Therefore 

^ < 2iVsh; < iVsh, (35) 

Lro Cjq 

Jia et alM computed Eq. (j34p for a Ag|YIG interface by first principles. The average of different 
crystal interfaces 

= 2.3x lOi^^^-im-2, (36) 

is quite close to the Sharvin conductance of silver (A^sh^o ~ 4.5 x lO^^fi"^ m^^). 
For comparison with experiment we have to include the Schep drift correction:— 

T 1 1 

(37) 



which leads to 

G,. f« 3.1 X lO^^f^-V-^ (38) 

One should note that the mixing conductance of the PtjYIG interface can then be estimated to be 
Gr ~ 10^^ m~^ since the Pt conduction electron density and Sharvin conductance are higher 
than those of noble metals. 

Using parameters p = = 8.6 x 10^'' Jim, = 12 nm, and A = 1.5 nm,'2S we see that the 
absorbed transverse spin currents with Gr = Gr and Gr = G™'^^ obtained from above for a Ag|YIG 
interface are 44% and 70% of the value for a perfect spin sink G,. — )• oo, respectively. For a Pt|YIG 
interface this value should be even larger. 

In order to compare our results with the observed SMR, we have to fill in or fit the parameters. 
The values of the spin-diffusion length and the spin Hall angle differ widely.— In Fig. [3] we plot the 
SMR for three fixed values of Gr- We observe that the experiments can be explained by a sensible 
set of transport parameters (G^, A, ^sh) that somewhat differ for the two representative samples 
reported in Ref. |25|. Generally, the SMR increases with a larger value of Gr but decreases when A 
is getting longer. These features are in agreement with the discussion of the simple limits above. 
Sample 1 in Ref. [2^ has a larger resistivity but a smaller SMR (ratio), implying a smaller spin 
Hall angle and/or smaller spin-diffusion length. When we fix the spin Hall angle 0sH = 0.06 and 
the spin-mixing conductance G,- = 5 x 10^^r2~^m~^, the corresponding estimated spin-diffusion 
lengths of Samples 1 and 2 are Ai ~ 1.5 nm and A2 ~ 3.5 nm, respectively. 

Finally we discuss the AHE equivalent or SHAHE. From experiments Ap2/p ~ 1-5 x 10~^ for 
p = 4.1 X 10"'''^ m and = 7 mnM Choosing = 0.05, A = 1.5 nm, and G^ = 5 x m'^, 
we would need a Gi = 6.2 x lO^^ri"^ m~^ to explain experiments, a number that is supported by 
first principle calculations 
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IV. SPIN VALVES 



In this section we discuss F(m) |N|F(m') spin valves fabricated from magnetic insulators with 
magnetization directions ifi and m' . The general angle dependence for independent rotations of m 
and m' is straightforward but tedious. We discuss in the following two representative configurations 
in which the two magnetizations are parallel and perpendicular to each other. We disregard in 
the following the effective field due to Gi such that the parallel and antiparallel configurations 
m = ibm' are equivalent. Moreover, we limit the discussion to the simple case of two identical F|N 
and N|F interfaces, i.e., the spin-mixing conductances at both interfaces are the same. 



A. Parallel Configuration (to • m' = ±1) 

When the magnetizations are aligned in parallel or antiparallel configuration, the boundary 
condition ji^\dN) = —ji^^ applies. We proceed as in Sec. Illll to obtain the spin accumulation 



y + rh X (m x y) 



2XGr tanh ^ 
a + 2XGr tanh 
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and the spin current 
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The spin currents at the bottom and top of N are absorbed as STTs and read 



.■SH 
JsO 



.•SH 



m X (rh X y) 



2XGr tanh ^ 



2A 



a + 2XGr tanh ^ 



leading to opposite STTs at the bottom (Tgtt ) and top {r^tt ) ferromagnets 



(39) 



(40) 



since jP{rh) = (0) = (c^at) = -j^s"" ' (rh'). 

The longitudinal and transverse (Hall) charge currents are 



^(B) 
'stt 



•(F) ^ 



2e 



Ji^^(O) 



' 'stt 



Jcjlong 
Jc, trans 



1+ 



'SH 



1-1 



m: 



2XGr tanh ^ 
^' a + 2XGr tanh ^ 



-Og-^rrixmy 



2XGr tanh ^ 



cosh 



a + 2AG^ tanh If cosh |f 



and the longitudinal and transverse resistivities read 

Plong = P + Apo + Api (l 

Ptrans = Apim^niy, 

where 



P 
Api 

P 



2 2A dN 
-^s^d^'""^2A' 
9^^ 4A2Grtanh2|f 
dN a + 2XGr tanh ^ ' 



cosh If 



(41) 



(42) 
(43) 

(44) 
(45) 

(46) 
(47) 
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FIG. 4: (Color online) Calculated Api/ (^6*1^) in an F|N|F spin valve as a function of spin-diffusion length 
A with (In = 12 nm, = 5 x 10^"* fi"^ m~^, and p = 8.6 x 10~''r2m chosen from Sample 1 in Ref. [25i 
Api/ (p^'sh) in an N|F bilayer is plotted as a dotted line for comparison. 



Figure m shows Api/ (p^gjj) with respect to the spin-diffusion length in an F|N|F spin valve with 
parallel magnetization configuration. Compared to N|F bilayers, the SMR in spin valves is larger 
and does not vanish in the limit of long spin-diffusion lengths. 



B. Limit X/cIn > 1 



The spin accumulation for weak spin-flip reads 
leading to the spin current 



(T -|- djSfGr 



rh X (m X y) 



2z — d]\f 



fs A/djv>l d^Gr 



^•SH 



-rh X [rh x y) . 



a -\- djsfGr 

In contrast to the bilayer, we find a finite SMR in this limit for spin valves: 



Jc.long A/djv>l 
Jc.trans A/q!jv>1 

IT 
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n2 
^SH 
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dN Gr 
a + djyGr 
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ml) 



rj2 
^SH 



djv Gr 

a + dNG. 



-m^my 



-Ol-^ruxmy 



1 + el^m 



or 



(48) 

(49) 

(50) 
(51) 

(52) 

^ CJ -|- djSfGr ^ ^ 

Here we find the maximum achievable SMR effects in metals with spin Hall angle ^SH by taking 
the limit of perfect spin current absorption. Clearly this requires spin valves with sufhciently thin 
spacer layers. We interpret these results in terms of spin angular momentum conservation: The 
finite SMR is achieved by using the ferromagnet as a spin sink that suppresses the back flow of 
spins and the ISHE. This process requires a source of angular momentum, which in bilayers can 
only be the lattice of the normal metal. Consequently, the SMR is suppressed in the F|N system 
when spin-flip is not allowed. In spin valves, however, the second ferromagnet layer can act as a 
spin current source, thereby allowing a finite SMR even in the absence of spin-flip scattering. 
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djSfGr Gr'>o/dM 
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C. Perpendicular Configuration (m • m' = 0) 

We may consider two in-plane magnetizations rh = (cos a, sin a, 0) and m' = (— sina, cos ce, 0), 
which are perpendicular to each other. When a = 0, the first layer maximally absorbs the SHE 
spin current, while m' is completely reflecting, just as the vacuum interface in the bilayer. For 
general a: 

2A0, ^^^^^^ 

Ai? a + 2XGr coth ^ \ sinh ^ sinh ^ / 

fisy{z) sinh^^^ 2XGr / cosh „ cosh f .„ \ 

' = :r- cos^ a ^ sm^ a , (55) 

/^^ sinh^ cj + 2AG^coth^ \^ sinh^^ sinh % J 

f^sziz) = 0, (56) 
which leads to the components of spin current normal to the interfaces 



2XGr tanh ^ ( sinh sinh f 



Qxi — J 1 J H -r- cos a sin a, (57) 

jfo + 2AG^ coth ^ \ sinh ^ sinh ^ / 

7,„(z) cosh - cosh ^ 2 AG,, tanh 4f /sinh^^ „ sinh f „ \ , ^ 

i.T cosh If a + 2AG^ coth ^ sinh^ sinh ^ J 

The total current is the sum of those from the two ferromagnets at the top and bottom; in contrast 
to the parallel m = itm' configuration, they do not feel each other. We can extend the discussion 
from the previous subsection: the second F can be a spin current source, and we can switch this 
source on by rotating the magnetization from perpendicular to (anti) parallel configuration. 
The longitudinal and transverse electric currents read 

Jc,ion,{z) _ , coshg^ 2AG.tanht / sinh ^ ^mhf , \ 

77; — L -\- Cqtt 1 h l/Qxi -, — -, COS CX — 1 — bUl Ct , 

Jc° coshff (j + 2AG,coth^ \ sinh^ sinh^^ / 

(59) 

7ctrans(^) 9 2AGr tanh /sinh^^ sinh f \ , , 



jj^ a + 2AGrCoth^ ^ sinh^if sinh^^f ^ 

Since the angle-dependent contributions vanish upon integration over z, there is no magnetoresis- 
tance in the perpendicular configuration. 



D. Controlling the spin-transfer torque 



Like the SMR, the STT at the N|F interface depends on the relative orientation between rn 
and m', too. We may pin m = x and observe how the STT at the bottom 
changes with rotating m' = xcosa + ^sina. Figure O displays the ratio /3 defined as 



/3(a) 



-(B) /-^ M 



rg\-t (x, X cos a + y sin a) 



'stt 



{x, x] 



(61) 



as a function of a for some spin-diffusion lengths. Only when X <^ d^, /3 remains constant under 
rotation of m'. A larger spin-mixing conductance and smaller d]\[ enhances the SMR as well as 
angle dependence of f3. This modification of the STT should lead to complex dynamics of the spin 
valve in the presence of an applied current and will be the subject of a subsequent study. 
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FIG. 5: (Color online) The ratio /3(a) which characterize how Tg\.^ changes with respect to the relative 
orientation between rh and rh' . We adopt the transport parameters = 12 nm, p — 8.6 x 10~''r2m, and 

V. SUMMARY 

We developed a theory for the SMR in N|F and F|N|F systems that takes into account the 
spin-orbit coupling in N as well as the spin-transfer at the N|F interface(s). In a N|F bilayer 
system, the SMR requires spin-flip in N and spin-transfer at the N|F interface. Our results explain 
the SMR measured in Ref. [i^ both qualitatively and quantitatively with transport parameters 
that are consistent with other experiments. The degrees of spin accumulation in N that can 
be controlled by the magnetization direction is found to be very significant. In the presence of 
an imaginary part of the spin-mixing conductance Gi we predicted a AHE-like signal (SHAHE). 
Such a signal was observed in Ref. [3Q and can be explained with values of Gi that agree with first 
principles calculations.^^ We furthermore analyzed F|N|F spin valves for parallel and perpendicular 
magnetization configurations. A maximal SMR ~ ^SH found for a collinear magnetization 
configuration in the limit that the spin-diffusion length is much larger than the thickness of the 
normal spacer. The SMR vanishes when rotating the two magnetizations into a fixed perpendicular 
constellation. The SMR torques under applied currents in N are expected to lead to magnetization 
dynamics of N|F and F|N|F structures. 
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